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Abstract 

Tuning interactions in the spin singlet and quintet channels of two col¬ 
liding atoms could change the symmetry of the one-dimensional spin-3/2 
fermionic systems of ultracold atoms while preserving the integrability. Here 
we find a novel SO (A) symmetry integrable point in the spin-3/2 Fermi gas 
and derive the exact solution of the model using the Bethe ansatz. In con¬ 
trast to the model with 577(4) and 50(5) symmetries, the present model with 
50(4) symmetry preserves spin singlet and quintet Cooper pairs in two sets 
of SU (2) <g) SU (2) spin subspaces. We obtain full phase diagrams, including 
the Fulde-Ferrcl-Larkin-Ovchinnikov like pair correlations, spin excitations 
and quantum criticality through the generalized Yang-Yang thermodynamic 
equations. In particular, various correlation functions are calculated by using 
finite-size corrections in the frame work of conformal field theory. Moreover, 
within the local density approximation, we further find that spin singlet and 
quintet pairs form subtle multiple shell structures in density profiles of the 
trapped gas. 
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1. Introduction 


10] and the two-orbital magnetism of SU(k ) sym- 
the SU(2) <g) SU(6) symmetry fermionic atoms of 


Large-spin atomic fermions displaying rich pairing structures and diverse 
many-body phenomena could be realized through controlling interactions in 
spin scattering channels. Promising theoretical progress has been made on 
large-spin atomic fermions in the context of multi-particle clustering super¬ 
fluidity Eli SB- In particular, fermionic alkaline-earth atoms can exhibit 
an exact SU(k ) symmetry with n = 2/ + 1 [s], 0], where / is pure nuclear spin. 
Experimental explorations of these fermionic systems have been reported on 
the trirner state of 6 Li atoms [[7, §j, 0 ], the /SC/(10) symmetry fermionic gas of 
87 Sr atoms with / = 9/2 
metry in 8 'Sr atoms [111 ]. 

173 Yb with its spin-1/2 isotope |l2|, the SU(6 ) Mott-insulator state of 173 Yb 
atoms EB arif l two-orbital Hubbard model EB> e ^ c - Despite much theoreti¬ 
cal and experimental efforts on the large spin atomic systems, understanding 
spin pairs and large spin magnetism is still rather elusive |l5|, see a recent 
review flhj ]- 

In contrast to the conventional spin-1/2 electronic magnetism [1 
one-dimensional (ID) large-spin ultracold atomic fermions exhibit richer high 
spin phenomena Es, 2ol [. Preliminary study of multi-colour superfluidity and 
quartetting ordering 21, 22| was also carried out through few simplest inte- 
grable large spin fermionic systems, such as the SU(A)- and S'0(5)-invariant 
spin-3/2 fermions [ 23 ], 24, 25]. However, there is still little progress toward 


E3 Es 


to the understanding of large spin pairing in non -SU(k) symmetry models 
26‘, 27, 28|. To this regard, the spin-3/2 fermionic system comprises an ideal 


model towards to the precise understanding of large spin lion-SU(n) symme¬ 
tries. From experimental point of review, spin-3/2 systems can be realized 
with alkali atoms of 6 Li, 132 Cs and alkaline-earth atoms of 9 Be, 135 Ba and 
137 Ba " 


19]. see the recent experiment 29] for more details 


Advances in controlling ultracold atoms provide us promising opportuni¬ 
ties of realizing quantum many-body systems with more tunable parameters 
30 . 3l| . For example, S-wave scattering in these spin-3/2 atomic fermionic 


systems acquires six effective interaction channels according to the Clebsch- 
Gordan coefficients. One is the spin singlet channel with total spin-0 and 
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the others are the spin quintet channels with total spin-2, while the S-wave 
scattering in the channels of total spin-1 and -3 is forbidden. When the in¬ 
teractions in all scattering channels are the same, the system exhibits SU (4) 


symmetry and is integrable in ID [23j, 24|, see model (i) in Table I. However, 


if the interaction in quintet channels is different from that in the singlet chan¬ 
nel, then the system can exhibit a 50 (5) symmetry without fine tuning fljj ]. 
see model (ii) in Table I. The model with 50(5) symmetry is integrable at a 
special point 


Table 1: Three integrable points for spin-3/2 quantum gases. 


No. 

9oo 

92,2 

92,1 

92,0 

92,-1 

92,-2 

i 

c 

c 

c 

c 

c 

c 

ii 

3c 

c 

c 

c 

c 

c 

iii 

c 

—c 

c 

—c 

c 

—c 


In fact, the scattering lengths in the quintet channels may be tuned to 
breakdown the SU( 4) or 50(5) symmetries while preserving the integrability. 
In this paper, we propose an integrable spin-3/2 atomic Fermi gas with 50(4) 
symmetry, see model (iii) in Table I. Based on the exact solution, we find that 
the model has subtle spin pairing phases and exhibits Fulde-Ferrel-Larkin- 
Ovchinnikov (FFLO) like pair correlations. For 5f/(4) symmetric model, 
there are no pairs when the interactions are repulsive, while for attractive 
interactions, four-body bound states appear in the ground state 23, 24 


With the competition between chemical /i and external magnetic held h , 
three-body bound states and two-body pairs also could emerge in the ground 
state 0- However, for the present 50(4) symmetry model with either 
repulsive interaction (c > 0) or attractive interaction (c < 0), bound pairs 
could emerge in the ground state regardless of the sign of c. In contrast to 
the 50(4) case, we also find that there are no four- and three-body bound 
states in the 50(4) system. We further show that the integrable 50(4) 
model presents universal quantum criticality with dynamic exponent z = 2 
and correlation exponent v — 1/2 in terms of different pairing states. 

The paper is organized as following. In Section [2j we discuss the sym¬ 
metry and conserved quantities of the integrable 50(4) model in the frame 
work of Yang-Baxter equation. The exact solution, thermodynamic limit and 
thermodynamic Bethe anatz equations are derived in Section [3J In Section 
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SI ground state properties and full phase diagrams are discussed in grand 
canonical ensemble. The elementary excitations of this model and Luttinger 
liquid behaviour are studied in detail in Section [5] In Section [6j various 
asymptotics of correlation functions for several quantum phases in strong 
coupling limit are calculated explicitly. The equation of states and the quan¬ 
tum criticality of the system are further studied in Section [71 followed by 
discussions in Section El 


2. Integrable S l O(4) model and conserved quantities 

We consider the model Hamiltonian 
n d 2 n 

^ = - Xj )- hM, (2,i) 

j =1 3 i^j Im 


that describes dilute spin-3/2 atomic gases of N fermions with contact in¬ 
teraction constrained by periodic boundary conditions to a line of length L. 
Here we let 2 m — H — 1. The last term in the Hamiltonian (12.Tj) is the Zee- 
man energy. The spin polarization (magnetization) is given by M — /?, 


where f z is hyperhne spin moment in he ^-direction. The projection opera¬ 
tor Pi™ = \lm)(lm\ projects the total spin-/ state onto the spin-m state in 
the ^-direction of two colliding atoms i and j. In the above equation, the 
summation YPm * s carried out for / = 0, 2 and m = —/,—/ + 1, • • • , Z. The 
interaction strength in the channel |Zm)(Zm| is given by gi m = —2Z?, 2 /(ma/^), 
where 


n lm - 
a i d ~ 


2 a lm 


Am 


1 -c- 


a _l 


is the effective ID scattering length depending on the 3D scattering length a lm 
in the channel |Zm)(Zm| [32|. The constant C ~ 1.4603 and (Z,m) stands for 
the interaction channel. Remarkably, the model (12 .1 [) exhibits three classes 
of mathematical symmetries that preserve the integrability by appropriate 
choices of the interaction strength via c lm = —2 /dfp = 7ngi m /h 2 , see Table 

m 

For model (iii), two sets of interacting channels {(2,1), (0, 0), (2, —1)} and 
{/2, 2), (2, 0), (2, —2)} form two spin subspaces with interaction potentials 
Hi = /§' + R 1 ;-" + Pff 1 and H 2 = P-f + P-f + P‘-f\ respectively. The full 
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interaction potential is 

Vii 

= 0^ 

-cV 2 , 
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also be written as 
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0 J 


Consequently, for c 

> 0 (or c < 

0), the potential V 2 (or 

Vi) presents attractive 


channels while V\ (or V 2 ) is repulsive. 

The multipole operators are defined by frii ] 

rytl _ \ ' rfil ijil—1 l,m—m' 

J-m / y 1 m' 1 m—m'^l, m > 

)7l'=0,±l 

m = l, l — 1, • • • , —l, l = 2, 3, 


i? = -(/* + i/')/V2, Tg = r, T-, = (/* - i/ ! ')/V2, (2.3) 

where = (ji, mi; j 2 , rri21J, m) are Clebsch-Gordan coefficients and 

/ x,y ’ 2 are the spin operators of the spin-f atoms 
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There are 15 multipole operators and only following 6 are commutative with 
the interaction potential 


/ 1 0 0 0 \ 


rp 3 _ 

i 0 “ 


2v'l0 


0 
0 

VO 

/ 0 0 


T\ = - 


T? = V 3 


-3 0 0 

0 3 0 

0 0 - 1 / 

0 0 \ 

0 0 


0 0 

10 0 0 
\ 0 -1 0 0 / 

/ 0 -1 0 0 \ 

0 0 0 0 

0 0 0 1 


\ n 


n 


rp3 _ 3 
2 “ 2 


T f 


1 

2 


T 2 , = 


/ 0 
0 
0 

V0 

/ 3 0 
0 1 
0 0 
V o 0 
/ 0 
1 
0 


0 1 0 \ 
0 0-1 
0 0 0 
0 0 0 / 

0 0 \ 

0 0 
-1 0 
0-3 / 


0 0 0 \ 

0 0 0 

0 0 0 

0-10/ 


.(2.5) 


Using above six multipole operators and the Pauli matrices, we define the 
pseudo spin-operators as 


J 1 '* = gbi + Tl) = iv <8 V, f’* = ^=(fy - f, 2 ) 

j'- y = jhi - h s ) = 1^ ® » z , f- y = + f?) 


d 2 <8) d®, 


—d 2 <g) dU 


i 1 ’* = \ 1 2f„ 2 - yy f„ s /. 


\/To r 

3 


1 

2 ( 


i V = l lTf + b^^/3] = ^®*‘ 


2\/lO r 
3 


1 

2 J 


( 2 . 6 ) 


The operators j d ' a with d — 1,2 and a = x, y, z satisfy the commutation 
relations 

= (27) 

7 

where e a/ g 7 is the fully antisymmetric tensor. We clearly see that j 1,a and 
j 2 ’" generate two SU( 2) Lie algebras and the six operators j d,a form the 
generators of the .SO (4) Lie algebra. 
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The attractive potentials may lead to spin-J pairs with J = 0, 2 in the 
quasi-momentum space. Consequently, the generators of these spin-J pairs 
comprise the following two sets of representations regarding to the attractive 
channels V 2 and V\ 


02,0 — -y= ('0-3/203/2 _ 01/20-1/2), 


02,2 — 01/203/2, 

02,-2 — 0-3/20-1/2, 

(2.8) 

00,0 — (0-3/203/2 + 01/20-1/2), 


o> 

to 

■€-> 

i 

co 

to 

02,-1 — 0-3/201/2- 

(2.9) 


Using the language of second quantization, the pseudo spin operators can 

be expressed as j d,a = xpj d ’ a, tp, where xp = ( 03 / 2 , 01 / 2 , 0 - 1 / 2 , 0 - 3/20 and 
the superscript t denote the transposition. One can check that following 
commutation relations are valid 


[i 1±,2 > 02, 1 ] = 

[j 2±,Z , 02 , 2] — 0 , 


[i 1±,2 >0 2 ,- 1 ] = °> [i 1±,2 >0o,o] = °> 
[] 2± ’ Z , 02,—2] = 0 , [ j 2± ' Z , 02,0] = 0 . 


( 2 . 10 ) 


The pseudo spin operators are commutative with the interaction potential. 
Therefore, the model has 5*0(4) symmetry. This symmetry does not preserve 
the number of atoms with each internal degree of freedom. But the total 
number of atoms are conserved. Meanwhile, the quantities 


J 3/2 — N 3/2 ~ JV-3/2, Jl /2 — N 1/2 ~ N- 1 / 2 , (2-11) 


are conserved. Here Ni is the number of atoms with spin-i component. The 
spin polarization along the ^-direction, which is given by M = | J 3 / 2 + | J 1 / 2 , 
is also conserved. 


3. The Bethe ansatz solution and thermodynamic limit 

The model (12. ljl with 50(4) s ym metry can be solved by means of the 
coordinate Bethe ansatz (BA) [331, ll7|. The BA wave function of the model 
reads 

'P(x 1 a 1 ,x 2 <T 2 , ■ ■ ■ ,x N a N ) = (3.1) 

V,Q 
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where Xj and aj are the coordinate and the spin in ^-direction of the j-th atom 
respectively. Here kj are pseudo-momenta of the particles, j = 1, 2, • • • , N, 
and Q and V are the permutations of { 1 , 2, • • • , N}. O(Q) is a series of pro¬ 
duction of step functions, i.e. O(Q) = 9(xq 2 — xq 1 )9(xq 3 — xq 2 ) ■ ■ ■ 9(xq n — 
xq n _ 3 ). As usual, 9(x) = 1 when x > 1 and 9(x) = 0 otherwise. In particular, 
the key ingredient of BA is the superposition coefficients of these plane waves 
A(Q, V), which can be consequently determined by the two-body scattering 
relation below. 

The wave function (13.111 of the many-body interacting fermions (12. ip 
should satisfy the fermionic statistics and time independent Schrodinger 
equation HE = E#. These restrictions lead to the scattering relation 
between two superposition coefficients A(Q( ba \ r p9 ,a '>) = SQ at Q b (k'p a — k Vb ) 
A(Q < - ab \'P < - ab ' > ), where the two-body scattering matrix 


sm = +ZY, kr, 

k — ic k + ic 

2=1,3 m=—l 

satisfies the Yang-Baxter equation 

S ab (X)S ac (X + u)S bc (u) = S bc (u)S ac ( A + u)S ab (\), 


(3.2) 


(3.3) 


which guarantees the integrability of the model. In the above equations 
we denoted Q (a6) = {■■■, Q a - i, Q a Qb , Qb+ i, ■ ■ ■ }, Q {ba) = {■■■, Qa-i, QbQa, 
Q b - |_i, • • • }, and V^ and V^ have a similar dehnition. Eq. (13. 2 p indicates 
that there is no interaction in the total spin-1 and -3 channels. This is the 
consequence of the symmetry of the wave function. The periodic boundary 
conditions give rise to the following eigenvalue problem 

eikjL Sj,j-i(kj - i) • • • S jt2 (kj - k 2 )S jtl (kj - k x ) 

xS j:N (kj -k N )--- S jJ+1 (kj - k j+ i)A = A, (3.4) 


which is used to determine the quasi-momenta k s. 

By using the nested algebraic BA 34], [35] and after complicated calcula¬ 
tions, we obtain the energy of the system as E — kj — hM, where the 













Figure 1: (Color online) Spin-flipping process: for c > 0, the spin rapidity A provides 
the spin flip with a total spin change of 1 whereas v gives the spin flip with a total spin 
change of 2. For c < 0, the rapidity v provides the spin flip with a total spin change of 1 
whereas A leads to a total spin change of 2. 


quasi-momenta {&;} should satisfy following BA equations, 

Mi M2 

e 1 kiL = n ei (hi - Xj) e_i (k t - z/,), i = 1 • • • N, 
3 = 1 i=i 

TV Mi 

J3 e ± (Aj - fci) = - JJd (Xj - Xe), j = 1 • --M u 

i= 1 t=\ 

N Mi 

13 e 1 (vj - fcj) = - JJ e 1 (i/j - u e ), j = 1 • • • M 2 . 

i = 1 " £=1 


(3.5) 


Here A and v are the spin rapidities, and we denoted the function e a (x) = 
x ~) ac . The BA equations coincide with the two-band model of electrons in 
ID [36|. From these BA equations we can obtain the full phase diagrams and 
thermodynamics of the spin-3/2 interacting fermions with 50(4) symmetry. 
In Eq. (13.51) . the quasi-momenta {/q} are nested with the spin rapidities A 
and v. For c > 0, the spin rapidity A provides the spin flip with a total spin 
change of 1 whereas v gives the spin flip with a total spin change of 2 , see 
Fig. [0 Thus the quantum numbers are given by M\ = N 1 / 2 + N_ 3 / 2 and 
M 2 = 1V_i/ 2 + iV_ 3/2 for {A} and {u}, respectively. For c < 0, the roles of 
the above v and A are swapped. Thus AA / 2 and iV_ 1/2 exchange in M\ and 
M 2 in this case. 

Finding root patterns of the BA equations (j3.5[) presents a big theoretical 
challenge towards to understanding the physics of the model. Here we find 
that the quasi-momenta {fcj} can be either real or complex conjugated pairs 
at the thermodynamic limit, L —» 00 and n = N/L is a constant. The 
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real roots k UjZ G M with z = 1, 2, • • • , N u characterize the quasi-momenta 
for unpaired fermions. The complex conjugated roots k p _ z e C with z = 
1,2,--- , N p denote the quasi-momenta of bound states, i.e., spin-J pairs 
with total spin J = 0, 2. We observe that the BA equations (I3.5[) present the 
spin pair bound states solely depending on the attractive interaction channels 
of Vj and V 2 (or say c < 0 or c > 0), see Eqs. (12.8p and (12.91) . Explicitly, for 
c > 0, the momenta for a bound pair of two-atom with different spins have a 
complex conjugate pattern k pz = k P:Z ±ic/2 with a binding energy e p = c 2 /2 
in the thermodynamic limit. Where a real spin rapidity A Pi2 = k PtZ is coupled 
together with k pz . So do the spin rapidities u p . z in the case of c < 0. 

For the ground state, we observe that the BA roots only involve real k, 
complex conjugated k, 1-string v and 2-string v for c > 0, see Section^ Here 
we denote the corresponding density distribution functions as p^\ p^ p \ p (l> 
and p^. Whereas, for c < 0, the 1-string and 2-string A populate the ground 
state instead of the u. At finite temperatures, the spin rapidities A’s and 
z/’s evolve different lengths of spin strings similar to the spin-1/2 Heisenberg 
integrablc model 


K -‘ + 2 

(n + 1 - 

2j)i|c|, 

Z= 1,2, •• • , Mi >n , 

1 

v n ,z + 2 

(n + 1 - 

2j) i|c|, 

z — 1, 2, ■ ■ • , M 2) „, 

t £ 

3 = 1,2, 

■■■n, 



(3,6) 


where M l n and ilf 2)n are the numbers of n-string of A and u, respectively. 
Thus we have Mi.n = Mi (/ = 1, 2). The BA equations in thermodynamic 
limit are obtained by substituting the string hypothesis into eqs. (13.51) . By 
taking the logarithm of the BA equations and using the relation 


x — 1C 

x + ic 


e iMH 


9 n (x) = 2arcta,n ( — ), 
\ncJ' 


(3.7) 
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when c > 0, we get the logarithmic form of the BA equations 


2vr I^\k uz ) 


27T&\k pz ) 


2rrl (1) (\ nz ) 


2 irl (2 \v nz ) 


Np OO M lrn 

k\izL y ^ 9 i ( k nz k p y) y ^ y ^ 9i!±(k xlz A my) 

y=l m= 1 y= 1 

oo M 2m 

+ ^ ^ y ^ 9™(k uz — Vmy), 

m= 1 j/=l 

N P N u 

2k pz L y ' ^py') ^ ^ (^p2 ^uy) 

y=l y=l 


OO M 2m 

+EE ^lm (^pz ^my ) 5 

m= 1 y=l 

AT U oo Mi m 

Efe^-u-EE- 4 nm(^nz 

y= 1 ra=l y=l 


^my) ? 


AT U 

^ ^ y) “1“ ^ ^ ^4.1n(^n2: ^py) 

y—i y=i 

oo Mi m 

^ ^ ^ ^ ^4-nra ( ^my ), (3,8) 

m= 1 y=l 


where ^4. nm — 0 n+m -f- 2 n Am —2 + ••• + 0 |n— m| —2 1 ~t~ 0 \n—m\ &n TTi^o* Horo the 

2 L 2 -2- J -2- ’ 

quantum number {/(^’(pMbd 2 ) j takes integers or half odd integers as the 
following 


fZ + - (iVp + M 2 ,in ), G Z + — (l + iV u + Ap), 

m m 

J (1) eZ+i(l + n + iV u ), J (2) G Z +t(l + n + AT). (3.9) 


The BA wave function acquires that any two rapidities of each branch are 
not equal. Otherwise the wave function is zero. According to this restriction, 
the quantum numbers Is of each branch take different numbers in parameter 
space [34!] • The solutions of the BA equations f!3.5p are classified by a set of 
quantum number |/( u )> (p)> (A ( 2 ) j. Thus the momentum of the BA eigenstate 
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is given by 


P 


2vr 

T 


E /(u> (*=».*)+E /<p, (^) 


+ ^/W(a„j + ^/‘ 2 >kj 

nz nz 


(3.10) 


Accordingly, the energy E and the magnetization M of the model are given 
by 


E = 

M = 


N u N p 

(Kz ) 2 + 2 + ^Vr - hM, 


3 =1 


3 =1 



+ 3N p — YXM ltn + 2M2 !H ), 

n 


(3.11) 

(3.12) 


respectively. 

At finite temperature, some quantum numbers Is (called vacancies) are 
occupied in parameter spaces whereas some of those quantum numbers are 
not occupied. The unoccupied BA roots are called unoccupied vacancies. In 
spin sector, the unoccupied vacancies are regarded as holes of the spin strings. 
The spin strings are nothing but spin wave bound states. For thermodynamic 
limit, i.e. L, N —> oo and N/L is finite, the strings and holes can be treated 
as density distribution functions p(k) and Ph(k), respectively. It follows that 
-A = p{k) + Ph(k). Consequently, we obtain the integral form of the BA 
equations in the thermodynamic limit, 


pi + p" = ^ - q * p" - £ * K’”‘ - p %m \ • 


Ph + P v = - hi * p“ - ai * p v + y A^ m * p*>"*, 

71 2 


£A 

m 

An | A.n y, . _u \ ^ Af . A,m 

Ph +P = «§ * P ~ 2_^ A nm * P > 


„2 ,m 


pT + p 2 ” = a f » p" + i' ln *p p -'E,A' r 


* n 2 ’ m 

nm * r 


(3.13) 


Here for c > 0, p 1,n and p 2,n are the string densities of n-strings of A and 71- 
strings of z/, respectively. Whereas for c < 0, p 1,n and p 2,n are the densities of 
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n-strings of v and n-string A, respectively. The * denote the convolution / * 
9(k) = f(k — k')g(k')dk'. The integral kernels are denoted as A' nm (k) = 

Anm(k) 5 nm (5(/c), A nm d\n—m\/2 T2(l(| n _ m |_|_2)/2 T ' ' 'T ‘Ad’(n+m—2)/2 T0(n+m)/2 

and a n (k ) = y • The p ln and p 1 ^ are the densities of particle and hole 

for the length ?7-strings in spin sector. The energy, momentum and total spin 
of the Hamiltonian are given by 

e = = j d kk 2 p u (k) + J &k(2k 2 — e p )p p (k) — hm, 

Jdkkp u {k) + J dk(2k)p p (k), 

Ad 3 ^^. . 

m = — = -n u + a 2 n p - Qi > y nm ln - a 2 2_^ nm 2n , (3.14) 


where o>i = 1, a 2 = 2 for c > 0, whereas = 2,a 2 = 1 for c < 0. We 
also denote 7i U)P = J d kp u,p (k) as the density of atoms in the scattering state 
and the density of the pairs, respectively. Moreover, m^ n = / d kp l,n (k) with 
/ = 1,2 are the densities of n-strings in the spin sector. 

At finite temperature T, following Yang and Yang’s grand canonical de¬ 
scription 38), the thermodynamic Bethe ansatz (TBA) equations can be 
obtained from minimization of the Gibbs free energy i? = E — TS — pN 
with S the entropy of the system. Here we define the dressed energies of the 
charge rapidities as E a (k) = Tin \p^{k) / p a (k)\ with a = u,p for unpaired 
fermions and pairs, respectively, thus the TBA equations are given by 


e u = k 2 - p n - ai/2 * £- + ^2 a n/2 * (e 1 ™ + £ 2 A), 

n 

e p = 2 k 2 — /i p — ai /2 * e" — a± * £ p + ^ A' ln * e 2 ”, 

n 

E ln = naih - a n/2 * e"_ - ^ A' nm * e 1 ™, 

m 

£ 2n = na 2 h + a n/2 * e" + A[ n * e v _ - ^ A' nm * e 2m , (3.15) 


where = — ln(l + p/ph ). The effective chemical potentials of unpaired 
fermions and pairs are defined by p u = p + (an + a 2 )h/ 2 and p p = c 2 /2 + 
2 p + a 2 h. For c > 0, e 1,n and £ 2,n are the dressed energies for n-strings of 
A and n-strings of n, respectively. Whereas for c < 0, e 1,n and £ 2,n are the 
dressed energies for n-strings of u and n-strings of A, respectively. 
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The TBA equations (I3.15P involve infinite branches of e l,n coupled to¬ 
gether. This imposes a big theoretical challenge to find the exact solutions 
of the TBA equations (13.151) . We observe that the dressed energy of the spin 
string tends to be a constant when the spin string length increases. The 
larger strings contribute the smaller energy as the temperature decreases. In 
order to discuss the solutions of the TBA equations, we present the TBA 


equations (13.15j) in the recursive form 

e u = k 2 — fi — G *£*?_ + G * + G * e 2 ! 1 , (3.16) 

£ p = 2 (k 2 - p- ^ - £ 2 ’\ (3.17) 

e 1,1 = — G * + G * (3.18) 

£ 2,1 = G * e" + G * e+ 2 , (3.19) 

£ 2,2 = G * £ p _ -\- G * + G * e+ 3 , (3.20) 

e 1>n = G* si” -1 + G * ^ n+1 , n> 2, (3.21) 

e 2 ’ n = G* 4 n_1 + G * ej n+1 , n > 3, (3.22) 

g.l,n £.2,71 

lim -— = ceih, lim 5— = a 2 h. (3.23) 

n—y oo ft n—»oo 77 , 


Here the convolution kernel is given by G(k) = (l/2|c|)sech(7rfc/c) and the 
function £ + is defined by £+(k) = Tln(l + p h /p). The thermal potential 
per unit length is p = p u + p p and the effective pressures are given by 
p u,p = — f^° oo dk£-(k) with r = 1 (r = 2) for the unpaired fermions (the 
bound pairs). The TBA equations (13 .1 61) - (13 .23 f) provide not only ground 
state properties at zero temperatures but also full thermodynamics at fi¬ 
nite temperatures. In the following Sections, we will solve the complicated 
TBA equations and discuss the physics of the spin-3/2 fermionic system with 
AO (4) symmetry. 

4. Quantum phase diagram and pairing signature 

In this section, we will study the phase diagram and the pairing signature 
of the model. A rigorous way of finding the ground state in grand canonical 
ensemble is to take the zero temperature limit in the recursive TBA eqs. 
(13. 161) - (13.231) . In the limit of T —> 0, we have £+{k) > 0 and £-{k) < 0. For 
a nonzero magnetic field h and c > 0, from Eqs. (13.231) . we observe that the 
large strings of A and v always have positive dressed energies. Therefore there 
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are no such strings in the ground state. From Eqs. fj3.2ip - fj3.22p . it is obvious 
that £ 1,n (n > 2) and e 2,n (n > 3) are positive so that the corresponding 
strings do not populate in the ground state. Even the magnetic held tends 
to zero, the dressed energies of these strings approach to 0 + . Therefore there 
are still no such strings in the ground state. The Eq. (13. 18p also gives a 
nonnegative dressed energy for the real A (or v) rapidity for c > 0 (or c < 0). 
Thus the dressed energy e 1,1 is also excluded from the ground state. 

In the T —> 0 limit, the TBA equations fl3.2ip - fl3.22p reduce to a new set 
of dressed energy equations that characterize the Fermi seas of the paired 
states and single atoms in terms of chemical potential p, and magnetic held 
h. The band filling of those Fermi seas with respect to p and h provide an 
analytical way to determine the full phase diagram in the / i-h plane. For 
c > 0, the dressed energy equations involve quasi-momenta k u and k p , and 
spin rapidities U\ and z/ 2 . The corresponding density distributions of them 
are denoted as p u ,p p ,p(b an d p( 2 ). For c < 0, the dressed energy equations 
involve k u , k p , zq and z/ 2 . We can define the ground state densities in a 
similar way as the ones for c > 0. The densities of these rapidities satisfy the 
following integral BA equations 


ph(k) + p(k) = p 0 - K*p(k), (4.1) 

where p h = (p£, p{, Ph\ PhYi P= (/A P p , P (1 \ P (2) )\ Po = (1/2^1/71,0,0)*. 
Here the superscript t means the transposition and integral kernel K is given 
by 


K 


( ° 

hi 

2 

—hi 

2 

\ -hi 


hi —hi 

2 2 

hi —hi 

—hi hi 

—hi — h 3 hi + h 3 
2 2 2 2 


-hi \ 


—a i — a 3 
2 2 

CL i (13 

2 2 

2ai + 5-2 



(4.2) 


The TBA eqs. fl3.16l) - fl3.23p in the T —)■ 0 limit reduce to the dressed energy 
equations which can be written in the following form 


e{k) = £q( k) — K * £~(k), (4.3) 

where e = (e u , £ p , ed), £ ( 2 ))t, £ _ — ( £ ^_ ) £ p _, £^\ e^)* and £ 0 = ( k 2 — pi, 2 k 2 — 
p 2 , cqh, 0 ( 2 ^)*. These equations are very convenient to analyze quantum 
phase diagram and quantum phase transitions in terms of chemical potential 
and magnetic held. 
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Figure 2: (Color online) Phase diagram of ID spin-3/2 Fermi gas: left panel for c > 0, 
right panel for c < 0. Contour plot shows the spin polarization rate M/N. The black solid 
lines divide the diagram into three quantum phases in quasi-monrentunr space: fully-paired 
phase (FP), fully-polarized phase (FF) and mixed phase of paired and unpaired fermions 
(FFLO). Here V stands the vacuum state. The red dashed lines divide the diagram into 
four mixed phases associated with different spin pairs, see in the text. Two sets of spin 
pair states in Eq. (12.81) give rise to the pairing phases in left and right panels, respectively. 


We observe that only the atomic pairs and the unpaired atoms populate 
the ground state. The absence of three- or four-body charge bound states 
is the consequence of the repulsive interacting channels due to the SO (A) 
symmetry. For c > 0, the existing pairs are </>\ ±2 and (j)\ 0 , whereas for c < 0, 
the paired states comprise (fr 2 ±1 and o- The competition among the pairing, 
chemical potential // and external magnetic field h gives rise to a rich phase 
diagram of this model. Following the method presented in J39|, full phase 
diagrams of the model can be both analytically and numerically worked out 
from the density equations (14. ip or from the dressed energy equations (14. 3 p . 
The phase diagram of the system with c > 0 is shown in the left panel in Fig. 
[2j For c < 0, different BA root patterns lead to a different phase diagram 
which is presented in the right panel in Fig. [2j 

In the following, we only consider the regime for c > 0. In Fig. [2j V 
stands for the vacuum. The phase FF denotes a fully-polarized phase of 
single 13/2) atoms in the region of strong magnetic field and small chemical 
potential. The FFLO-likc phase composes of the largest spin-2 component 
pairs 02,2 and excess fermions of 13/2) atoms that form two mismatched Fermi 
surfaces in quasi-momentum space. In this phase, atoms can be in either 
the scattering state 13/2) or bounded pair state 0 2 ,2- We show that spatial 
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Figure 3: (Color online) Density profiles of a spin-imbalanced ID Fermi gas of spin-3/2 
atoms in an harmonic trap with the setting (h/e p , N/y/e£) = (0,2.82), (0.6,1.75) and 
(1.8,3.32) for (a), (b) and (c), respectively, (a) The MP phase of ^>2,2, 02,0 and <j> 2,-2 pairs 
extends to a whole ID tub. (b) A core of the mixed phase III is surrounded by 6 wing 
shells of mixed states I, FFLO state and the pure 02,2 paired states FP. (c) A core of the 
mixed phase III is surrounded by 6 wing shells of mixed states II, the FFLO state and the 
fully-polarised phase FF, also see Fig. [2] 


oscillation of the pair correlation function solely depends on the mismatch 
between the Fermi surfaces of 13/2) and 11/2) atoms, i.e. A k F = 7r(n 3 / 2 — 
77,1/2). In the strong coupling region, the slowest decaying term of the pair 
correlation function G p (x,t) = {G\4>\ 2 ( a: j 2 ( 0 , 0)|G) is obtained, G p ~ 

Aq cos('K/S.kFx )/(\x + w u f| e,u |a; + iu p f| 0p ), where 9 U = 1/2, 6 P = 1/2 + n p /c 
and n UiP are the densities of unpaired 13/2) atoms and atomic pairs 02 , 2 , 
respectively. More detailed analysis will be further discussed in Section |6j 
Thus the momentum pair distribution has peaks at the mismatch of the 
Fermi surfaces A kp, which presents a characteristic of the FFLO phase. 

Moreover, the FP phase is the fully-paired state 0 2 , 2 > where the pair 
correlation function decays as a power of distance. The leading term of pair 
correlation function G p is given by G p = A/(\x + iv p t\ ep ) with 9 P = (1/2) (1 + 
I/ 7 ). This indicates that the pair correlation in the fully-paired phase has a 
longer correlation length than the one in the phase FFLO. According to the 
Luttinger liquid theory, the long distance behavior of the pair correlation of 
this type gives rise to the algebraic form G p ~ \x\~ 1 ^ 2K \ where the parameter 
K ss 1 — I/ 7 . However, the correlation function for the one particle Green’s 
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function decays exponentially. When the interaction is very strong, i.e. near 
the phase transition line from vacuum into the FP phase, the system behaves 


like a Fermionic super Tonks-Girardeau gas 40, 41 


The MP phase denotes a fully paired phase, but not a fully-polarized 02,2 
state. For c > 0, three kinds of spin pairs 02,2, 02,o and 02,-2 are involved 
in the MP phase. By increasing the chemical potential p and magnetic held 

h , the system may enter into three other mixed phases denoted as V, FF 
and FFLO Phases, see Fig. [21 The phase transition from FFLO phase into 
V phase occurs as the 2-string v emerges in the ground state. While the 
one from MP phase into V phase is driven by the appearance of real fc, 

i. e. involvement of pair breaking. The one from V phase into FFLO phase 
is driven by spin hipping processes involving the real v and the one from 
FF phase into FFLO phase is driven by spin hipping process involving the 
2 -string v. 

In a harmonic trap, quantum phase segments in density prohles can be 
used to identify different quantum phases in experiments with cold atoms. 
We can extract the threshold values of the phase boundaries in Fig. [2]through 
the density prohles of the trapped gas. Within the local density approxima¬ 
tion, the local chemical potential is replaced by /.i(x ) = p 0 — f° r the 

trapped gas. The total density of fermions is given by n = n u + 2n p . For 
hxed particle numbers, we plot the density prohles of the trapped gas in Fig. 
[3j In contrast to spin-1/2 Fermi gas [421,0,0, 45], spin quintet pairs lead 
to multiple shell structures in the density prohles. The trapping center can 
be a mixture of three types of spin pairs as well as unpaired single atoms 
accompanied by the different shells involving the fully-polarized spin-2 pairs, 
the FFLO state, and the FF phase, see Fig. [3] 


5. Elementary spin and charge excitations 

We consider the elementary excitations in the fully paired MP phase. 
Suppose that the external magnetic held is zero. The excitations in the charge 
sector of the system are similar to the ones in the SU (4) interacting fermions. 
However, the spin excitations in the spin sector are quite different from those 
of £77(4) Fermi gas. In the fully paired MP phase, the elementary charge 
excitation has an energy gap. In contrast to the SU (2) attractive Fermi gas 
0 |, there exist gapless low-lying spin excitations triggered by changing the 
quantum number of 2-string v in this paired MP phase. The change of the 
number of 2-string v results in necessary spin configuration changes in the 
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Figure 4: (Color online) (a) Dressed energies of different rapidities in the fully paired 
phase for h = 0 and 7 = 1 . The red solid line stands for the dressed energy of the real k 
(unpaired fermions). The black solid line denotes the dressed energy of the pairs (2-string 
of k). The blue dashed line shows the dressed energy of the 2-strings of v. The black 
dotted line is the dressed energies of higher strings of v. (b) Elementary excitation spectra 
for the model with 7 = 1 and h = 0 for the phase MP. One 2-string excitation of v splits 
to four spin -1 holes of v, i.e. creation of four spin -1 spinons. 


paired states of 0 2)2 , 02,0 and 02 ,- 2 - In particularly, breaking down a 2-string 
v leads to the total spin change of 4, i.e. creation of 4 spinons, see Fig. [5j 
Each spinon carries spin-1, and we denote this type of excitations as large 
spinon excitations. There is no such kind of spin excitation in the attractive 
5(7(2) Fermi gas. These spinons can be directly calculated from the BA 
equation (14.11) . One 2-string of v is excited that leads to four spin-1 holes of 
v. SM 2 = 5 /4. If one adds d holes of 2-string u into the ground state, the 
excited energy and momentum are given by 

A £ =-^y 2) (^ h ), AP = -^2 tt/' p (2) (z/| dz/, (5.1) 

i=i i =1 Jo 

respectively. Where z// are the positions of the added holes. For h = 0 and 
7 = c/n = 1, the 4-hole exciting spectra is shown in Fig. H(b). However, for 
c < 0, one 2-string A excitation splits into four spin-1/2 holes which give an 
excitation with total spin-2. 

I 11 order for our convenience to discuss correlation function, we will focus 
on three phases FF, FFLO and FP for the strong coupling regime and a finite 
external magnetic field h. These phases near the vacuum phase boundary in 
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Figure 5: (Color online) Three types of elementary excitations: particle-hole excitation, 
adding particles at different Fermi points and backscattering process. kp is the pseudo 
Fermi points at the right/left Fermi points, a = u, P denote the unpaired atoms and 
atomic pairs, respectively. The upper, middle and lower panels show the Type I, Type II 
and Type III, respectively. 


Fig. [21 It is important to note that in these three phases both spin rapidities 
of v and A are absent in the ground state. Thus the dressed energy equations 
(14.31) reduce to the similar structure as the one for the attractive spin-1/2 


Fermi gas 39, [47], 125]. However, spin excitations and finite temperature 


behaviour of these states in the 50(4) symmetry model are quite different 
from that of the attractive spin-1/2 Fermi gas. In general, at zero temper¬ 
ature, the Fermi points Q UjP are determined by the condition £ U (±Q U ) = 0 
and £ P (±Q P ) = 0, where Q U)P > 0. At the ground state, all the quantum 
numbers are filled up to the Fermi points, i.e. / u (fc u , 2 ) — J u (fc Ui2 _i) = 1 and 
I p (k P ,z) ~ I p {k PtZ -i) = 1 with k UyZ > k U)Z -\ and k PtZ > . All excitations 

can be classified as three types of excitations, i.e. particle-hole excitations, 
backwards scatting process and adding particles at different Fermi points. 
We denote these excitations as Type 1, Type II and Type III, respectively, 
see Fig. El 

The Type I elementary excitations are characterized by moving atoms 
(pairs) close the left or right pseudo Fermi point outside the Fermi sea, see 
Fig. El (1)- For example, we consider the excitations close to the right Fermi 
point Q u in real quasi-momentum space k. The process is to create IV+ holes 
at positions I u (k U:Nu ), I u (k u , Nu ) - 1, - • ■ , I u (k UtNu ) — N+ + 1 and to add N+ 
real ks outside the Fermi sea where the corresponding quantum numbers are 
I u (ku,N n ) + 1, I a (k U:Nu ) + 2, • • • , / u (fc u , 7 vJ + N+. Based on this setting, one 
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can perform standard calculation of finite-size corrections to the energy and 


momentum 


43,125 


This type of excitations can take place close to the 
left pseudo Fermi point — Q u with adding N~ holes below the Fermi point. 
In the branch of atomic pairs, particle-hole excitations also occur close to the 
left and right pseudo Fermi points with adding the quantum numbers N~ 
and N+, respectively. 

The Type II excitations arise from the changes of total number of unpaired 
fermions or bound pairs. This type of excitations are characterized by adding 
(or removing) atoms/pairs close to the Fermi points, see Fig. [5] (2). We 
denote the change of particle number as A1V U (or AN P ). 

The Type III excitations are created by moving particles from the left 
Fermi point to the right Fermi point and vice versa. This process is also 
known as backscattering, see Fig. [5] (3). These two types of excitations 
change the pseudo Fermi points. The quasi-momenta of single atoms/pairs 
with positive real parts are viewed as the right-going atoms/pairs whereas 
the ones with negative real parts are denoted as the left-going atoms/pairs. 
Both the Type II and the Type III excitations lead to the particle number 
difference between the right- and left- going atoms/pairs. We denote these 
number differences as 2 A D u and 2 A D p for the unpaired atoms and the pairs, 
respectively. 

Based on the above descriptions, all the three types of excitations can 
be unified in the following form of the finite-size corrections for the total 
momentum and energy of elementary excitations 


A P = 

2 vr 

T 

A E = 

2 t r 

T 


+] 


2 t r 


T 


) rrr 

— [AJV„A D a + N a AD a + A), - JV“], 


(5.2) 


E 

a(3 7 


-AN a [Z 1 ]p a vp[Z 1 ]^ 7 AiV 7 + AD a Z a pVpZ y pAD, 


E t '«( jv « + + jv .) 

a 

E" 


Z~'AN]l + (A DZ\\ + iV+ + JV-), 


(5.3) 


where the indices a, (3 ,7 = u, p and vp is the Fermi velocities of the cor¬ 
responding states. In the above equations Z = Z {w ^ is the dressed charge 
which is determined by 


z [w \k) = i (w) - k (w) * z (w \k ), 


(5.4) 
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and w stands for different phases, i.e. FF, FFLO and FP. For the fully- 
polarised phase, u>=FF and we have 

Z^ FF \k) = Z uu (k), J (FF) = 1, K (FF) = 0, Z^ = Z uu (Qu)- (5.5) 


For the pure paired phase, w=FP and we have 
Z^ FP \k) = Z pp (k), I (FP) = 1, K^ FP \k) = ai (k), Z( pp ) = Z pp (Q p ). (5.6) 
For the phase, ta=FFLO and we have 


^•(FFLO ) (k) = ( 

y Z pu (k) Z pp (k) 

T-(FFLO) ( 1 0 \ t-(FFLO) 

y 0 1 ) ’ 


K (FFLO) (fc) 


0 

ai(k) 


Z uu (Qu) Z up (Q u ) \ 
^pu(Qp) z pp (Q p ) ) 


ai(k) 

ai(k) 


(5.7) 


The value of dressed charges can be obtained by solving the dressed charge 
equation (15.4j) numerically. The asymptotic forms of the dressed charges in 
the strong coupling regime is given in Section [HI 

On the other hand, the excited energy and excited momentum can also 
be expressed by the conformal dimensions A± according to the conformal 


held theory [491, (50 


o 7 r 

A£ = t 5>„( AJ + A-), 

a 

ad = t E (A t - A °)+ 2 E AD ° k *, 


where kp )CC is the Fermi momenta, k-p a = 7r N a /L. It follows that 

a: + A; = N* + N- + jlZ-'AJV]^ + [ADZE, 
AJ - A; = JV+ - 1 V- + AN a AD a . 


(5.8) 

(5.9) 


(5.10) 


These relations provide us a analytical way to calculate the conformal di¬ 
mensions via fini te-size corrections 

2 A± = 2 N± ± AN a AD a + [A DZ] 2 a + ^[Z^AN}1 (5.11) 
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6. Asymptotics of correlation functions 

In this section, we will calculate various correlation functions of the sys¬ 
tem with strong interactions in three phases FF, FFLO and FP. At zero 
temperature, all correlation functions reveal a power law decay in the phases 
FF and FFLO. In the pure paired FP phase the pair correlation decay as 
a power law of distance whereas the single particle Green’s function decays 
exponentially. However, at finite temperatures, all correlation functions ex¬ 
ponentially decay in long distance. 

From the conformal held theory, the two-point correlation function for 
primary fields with the conformal dimensions A ± is given by j49| 

— A p — 27ri(52 c( N a AD a )x/L 

Go(x,t) = Y,Wl - ■ TT (6- 1 ) 

1 la 1A — i v a t) a [x + w a ty* a 

where Go(x,t ) = {G\0\x, f)O(0, 0)|G) is the correlator for the held op¬ 
erators 0(0,0) and 0^(x, t). Here the conformal dimensions A 1 * 1 are de¬ 
termined by (15.111) in terms of N ± , AN, AD and the dressed charge Z. 
For strong interaction regime, the leading order of the dressed charge is 
given by Z a p = S at p, where S at p is the Kronecker delta function. From this 
leading order of the dressed charge, we obtain the conformal dimensions as 
2A± = 2N± + (AD a ±AN a /2) 2 + 0(c~ 1 ). For the strong coupling regime, the 
order of 1/c corrections to the correlations are presented via the conformal 
dimensions 

2A± = 2 N± + (A D a ± ANJ2) 2 + + 0(c" 2 ), (6.2) 

where S a with a = FF, FP and FFLO are the hrst order corrections, see 
Table [2] 

In the fully-polarized phase FF, fermions are unpaired and there does 
not exist contact interaction. Thus the dressed charge reads Z" = 1. In the 
phase FFLO, the dressed charges are given by 



For the fully-paired phase FP we have the dressed charge of pairs Z£ = 
1-2Q p /ttc + (9(c- 3 ). 
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Based on the above settings, we may calculate the asymptotics of the 
correlation functions {G\0\x, t)O(0, 0)|Cr). For a given correlation function, 
the operator 0(x, t ) act on the ground state that gives the quantum numbers 
AiV Uj p. Expanding the operator 0(x,t ) with respect to the primary fields 
with conformal dimensions A± and their descendent fields, we may calculate 
the asymptotics of the correlation functions by the formula (16. ip . We present 
the conformal dimensions of various correlators in Table [21 In general, the 
real part of a correlation can be expressed as 

g(x) = Acos(2-kx /\ s ) \x — iv a t \~ 6a , (6.4) 

Q 

where X s = (^2 a N a AD a /L) _1 is the wave length of the oscillation and the 
exponent 9 a = 2A+±2A“. The Fourier transform of the correlation function 
near the Fermi velocity ko is given by [5l| 

g(k) oc [sign(A; - k 0 )} 2s \k - k 0 \ Vs , (6.5) 


where u s = 2A + + 2 A - — 1, s = — A” is ^ ie conformal spin and 2 s 

is always integer. 

The various correlation functions, for example, single particle Green’s 
function, pair-pair correlation, density-density correlation are listed in Ta¬ 
ble [21 The fully polarized case FF is regarded as a single component free 
fermionic gas. The correlation function of the field operator is given by 

G u (x,t) = (G'|^ /2 (x,f)'0 3 / 2 (O,O)|G). (6.6) 

This excitation associates with the quantum number A N u = 1 and A D u e 
Z+l/2. A few leading orders of the asymptotics of the single particle Green’s 
function are obtained by choices of the values of the AD, see the rows 1-3 in 
the Table [21 The leading term in the correlation function G u is given by 


G 


FF 


A, 


cos(7rn u x) 

n • 7T i 
\x — m u r| 


(6.7) 


with the setting A D u = ±1/2 and N+ = N~ = 0, see the row-1 and the 
row-2 in Table [21 Here n u is the density of single atoms and A 0 is some 
constant. The row-1 shows a left-going wave with A D u = —1/2, or say that 
adding an unpaired atom near the left Fermi point. The row-2 indicates a 
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Table 2: Various correlations functions for quantum phases FF, FFLO and FP. 
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right-going excitation wave with A D u = 1/2. We also consider the charge 
density correlation function G n (x,t ) = {n(x, i)n(0, 0)) in this phase, namely 


G n = (G\n(x, t)h(0, 0)|G) = n 2 + 


A nl cos( 27 t n u x) 


\x 


iv n t\ : 


+ 


A, 


n2 


\X 


iv u 


( 6 . 8 ) 


where A nl and A n2 are constants, see the rows 4-8 in the Table [2] 

In phase FFLO, both the spin-2 pairs and single spin-3/2 atoms coexist. 
The single particle Green’s function G u acquires A N u = 1, A1V P = 0 and 
AB UiP £ Z + 1/2. The few leading orders of the asymptotics of G u are 
indicated in the rows 9-12 in Tabled The first two leading terms read 

g fflo _ A i cos[7r(n u + n p )x] A 2 cos[7r(n u - n p )x] 

u \x + in u t | 0ul \x + iv^p 1 \x + in u t | e ' u2 \x + iv p t\ e P 2 ’ 

where Ai and A 2 are constants and the exponents read 


^1 = 1- 2 n p /c, 6 pl = 1/2 - 4n u /c - n p /c, 

0 u2 = 1 + 2?i p /c, 0 p2 = 1/2 + 4n u /c - n p /c. (6.10) 


In the correlation function 06.911 . the hrst term presents the result listed in 
the rows 9 and 10 in Tabled and the second term is presented in the rows 11 
and 12. For finitely strong interaction, the single particle correlation function 
decays slower than that of the free Fermions due to the correlations between 
pairs and between pairs and unpaired atoms. 

In the FFLO phase, the asymptotic pair-pair correlation is given by 

G p (x,t) = (G|0j)(£,f)</> P (O,O)|G), (6.11) 


with the quantum numbers A N u = 0, A1V P = 1, A D u £ Z + 1/2 and 
A D p £ Z. The leading order of the pair-pair correlation function is given 
explicitly by 


c _ H 3 cos[7r(n 3/ 2 - n l/2 )x] 
p \x + iv u t\ 9u3 \x + iv p t\ 9 P 3 ’ 

(6.12) 

where A 3 is a constant and 


0u3 = 1/2, 0 P 3 — 1/2 + n p f c, 

(6.13) 


see the rows 13 and 14 in Tabled Here 7i 3 / 2 and rti/ 2 are the densities of spin- 
3/2 and -1/2 atoms, respectively. We see that the long distance asymptotics 
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for the pair correlation function oscillates with the wave number 71 ( 713/2 — 
ni/ 2 ), a mismatch between the two Fermi surfaces of the two species. This is 
a characteristic of the FFLO-likc correlation. We would also like to address 
that the FFLO-like oscillation terms arise from Type 111 excitations, i.e. 
backscattering for bound pairs and unpaired fermions. The charge density 
correlation function G n (x,t) is presented in the rows 17-19 in Table [21 

In the phase FP, the pair-pair correlation function G p acquires the quan¬ 
tum number AN P = 1. We present the result of the long distance asymptotics 
of pair-pair correlations in the rows 24-27 in Table [2j The leading term of 
the pair-pair correction reads 


G « 

Orp ~ 


A d 


\x + iv p t\ e P 4 ’ 


a _ 1 , n P 

Opt-2 + - 


(6.14) 


where A 4 is a constant. This shows that the exponent 9 p 4 is smaller than 
that in the phase FFLO. In this phase the pair-pair correlation dominates 
the quantum correlation. Finally, we present the charge density-density cor¬ 
relations 


G n = {G\n(x,t)n( 0 , 0 )\G) = n 2 + A n 2 g(x,t), (6.15) 


with the correlation 


g(x,t) 


cos(27 rn p x) 
\x + iv p t\ dl ’ 


2 - 4n p /c, 


(6.16) 


where A n2 is constant. The related quantum numbers are listed in the rows 
28-31 in Table El 


7. Equation of state and quantum criticality 


The ID fermionic systems usually exhibit a rich resource of Luttinger 
liquids and show a novel magnetism in term of spin-change separation sce¬ 


nario, see recent review 46]. In order to understand large spin cold atoms 
in ID, it is very important to investigate the low temperature behavior of 
Luttinger liquids. We will prove that the Luttinger liquids of different pair¬ 
ing states comprise a universal low energy physics of the model. In general, 
the specific heat can be written in terms of sound velocities of the liquid 
phases, i.e. cl = (nTL/3)'^2 a l/v a , where the summation carries out over 
all the Fermi velocities in the charge and spin sectors. This result is valid for 
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the temperature below a crossover temperature T*, i.e. T -C T*. Here the 


crossover temperature can be determined from the equation of states 52 


Near a critical point and for the temperature T T*, the system lies in the 
critical regime, where the crossover temperature characterizes the energy gap 
T* ~ /i — /i c and/or T* ~ H — H c . In the critical regime, universal scaling 
behaviour of thermodynamical properties is expected to distinguish quantum 
criticality with the critical dynamic exponent z — 2 from the Luttinger liquid 
criticality of z — 1 [39]. 

For strong coupling regimes, only charge rapidities are left at the ground 
states for the phases FF, FFLO and FP. For low temperatures, the spin 
wave contributions to the dressed energies of these phases can be analytically 
derived. In the same fashion as the SU( 2) Fermi gas [39], we can express the 
TBA equations of the dressed energies as 


e a (k) = r a k 2 - * £_(*:) + f a (k), a = u.p. 


(7.1) 


Here r a is viewed as the effective mass of the corresponding charges. For 
the unpaired atoms, r u = 1 and for the pairs, r p = 2. In the above equa¬ 
tions, f a denotes the contributions of spin wave bound states to the dressed 
energies of the charge degree of freedoms. However, for h T, the spin 
contribution term f a will exponentially decay as e~ h ^ T and it is negligible 
in low energy physics. Follow the method proposed in 39], for the strong 
interaction regime, we find the following equation of states 


VaT 2 /A a 
Pa = - —Fi | — 


27T 2 


2 V T 


1 - 


D 


( 3 ) 


:/#/3 


c 3 r a rp 


)+0(c- 5 ). 


(7.2) 


Here the effective chemical potentials A a is given by 


A a Ha 


2/}W 

fi 


r p c 


PP 
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*i(£). 


of - 3 ^ v T 

ZTpTrze 6 


(7.3) 


and D ^ are constants resulted from the strong coupling expansion of the 
integral kernels K. Fj(A a /T) are Fermi-Dirac integrals which can also be 
written as polylogarithm functions 39]. D fd = F)( 3 ) = 0 in the phase FF. In 
the FLLO phase, we have 


= 


L) (3) = 


0 8 
8 1 


(7.4) 
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Figure 6: (Color online) Universal scaling behaviour of the density n at the quantum 
criticality for h/c: 2 = 0.55: (a) at the phase transition V-FF and (b) at the phase transition 


FF-FFLO. 


In the FP phase, these constants are given by D ^ = —/T 3) = 1. By iteration 
with p a and the effective chemical potentials /l UjP , we obtain a close form of 
the equation of sates for the FFLO phase as 


V 


Ff 1 + 


rjl — 

1 2 

~i=t i 

V / 27t 2 

rpi 

1 2 

+ —=Fi u 

2a/F 2 


1 + 


4 rf 
a/27t c 3 
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L 2 16 L 
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F ? 


.3 o 


+ 0(c 


-5a 


(7.5) 


where iy" = Fj(A a /T) and renormalized chemical potentials 7l UiP are given 
by 
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T 5/2 
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2-y/7 tc 3 
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X 5 / 2 
4 a/ 27 tc 3 


A /2 .(7.6) 


This spin-3/2 Fermi gas exhibiting rich quantum phase transitions pro¬ 
vides an ideal model to investigate quantum criticality of large spin inter¬ 
acting fermions. The phase diagram Fig. [2] has multifold critical points. In 
particular, the chemical potential and external held could drive the system 


29 


























from one phase to another as these parameter pass across the phase bound¬ 
aries in the phase diagram Fig. [2] at T = 0. All phase transitions in Fig. [2] 
are of the second order. Near the critical point, thermodynamical proper¬ 
ties evolve into certain universal scaling forms as temperature tends to zero. 
Below a crossover temperature (see the upper panel of Fig.[8]below) the quan¬ 
tum criticality of the Luttinger liquid gives the dynamical exponent z — 1. 
However, beyond the crossover temperature, a non-relativistic quadratic dis¬ 
persion leads to the free fermion criticality with z = 2, where new band 


excitations are involved 24, 25 


53 


For 7 > 1 and h T |/x 
the universal scaling function 


He |, the density and compressibility satisfy 


n - n0ti = r^+ 1 - 1 /-J)[(/r - 

k* - K* 0ii = T d / z + 1 - 2 / vz g i [{ l i - fi c i)/T l / vz l (7.7) 


where we find d/z + 1 — 1/vz = 1/2, 1 /vz = 1 with d = 1, critical exponent 
z = 2 and correlation length exponent v = 1/2. In the above equation, 
i — 1, 2, 3,4 denote the four phase transition boundaries of V-FF, V-FP, FF- 
FFLO and FP-FFLO. For i = 1,2, both the background density n 0:i and the 
compressibility /c 0 ,i are zero. For i — 3,4, we have 


^o,3 =--«> ^o,4 = -@(l - -P + \t3 2 \ 

2n 71 \ 71 7l z / 

1 1 2 1 / 3 6 2 \ 

k 0,3 — X- 5 ^0,4 o ( 1- P 5- 2@ )’ 

111 a 71 p \ 71 7l z / 


a = < 2 aAh — h P ) 1-|- (2ot\(h — hp)) 2 

l L 37 t 

( r 2 i 

(5 = < 2ai(h P — h) 1-1—(2 aAhp — h ))a 

l L 71 


(7.8) 


where hp = o>2e p /2 is the magnetic held at the four phases (V, FF, FP and 
FFLO) coexistence point. In the above equations the scaling functions Fi(x) 
and Gi(x) are 


J)(x) 

1 20F ’ 
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J):(a;) 

F_ 1/2 (2 x) 
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Figure 7: (Color online) Universal scaling behaviour of compressibility n at the quantum 
criticality for h = 0.05: (a) at the phase transition V-FP and (b) at the phase transition 


FP-FFLO. 


where the coefficients 


$1 = $ 2 = $i = $2 = 1 , 

-. , 2 1 2 ^ i 8 o l 7 02 

$3 — 1 a H — -a , $4 — 1 / 3 H — y/3 , 


9 7 ■tt / 1 o r 

7T 7T Z 7T 7T Z 


1 1 9 

$ 3 = H-a--" 2 

2vr 2vr 2 


10 


a 2 , $ 4 = 1 + -p - (7.10) 

7 r 7 t- 


Let us dehne the dimensionless density n = 7 i 0 T 1 / 2 and compressibility 
k* = then we have the scaling forms 

Ah = h - h 0>i = Jy(A n/T), A k* = k* - k 0ji = ^(A n/T), (7.11) 

where A/d = /i — /i ci . We observe that these physical quantities intersect at 
the critical point. This intersection nature can be used to map out quantum 


criticality of the model through the trapped gas at finite temperatures 39 


also see experimental study of the quantum criticality for the 2D Bose atomic 
gases [541, 55, 56, 57. 

We plot the critical properties of the density for the phase transitions 
V-FF and FF-FFLO in Fig. EE Whereas the critical properties of compress¬ 
ibility for the phase transitions V-FP and FP-FFLO are presented in Fig. [71 
Universal scaling behaviours of the density and compressibility read off the 
dynamical critical exponent z = 2 and correlation critical exponent v = 1/2. 
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Figure 8: (Color online) Upper panel: quantum criticality of the phase transition V- 
FFLO near the critical point fi c = — 1.5e p for the fixed h = e p and c > 0, where C denote 
the classical region, QC means the quantum critical regime and TLL is the Luttinger 
liquid. The white dashed lines indicate the crossover temperature. Lower left (right) 
panel: density (compressibility) vs ^ at different temperatures. The solid lines denote the 
analytical result from (17.1211 and the dot lines show the numerical solution from the TBA 
equations (18.151) . 


In Fig. H]we further demonstrate the quantum criticality near the phase 
transition from vacuum into FFLO at the fixed magnetic held h = 6b- This 
phase transition involves a sudden change of the densities of states of the 
spin-2 pairs and single spin-3/2 atoms. At low temperatures, we obtain the 
scaling functions 
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N, 


20F _ 
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F_ i 
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20F _ 


F 3 

2 


h he 
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h he 
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+ 22 F_ 


+ 2 2 F_ 3 
2 


2(h he) 

T 

2(/i — /T c ) 


(7.12) 


where fi c = — 1.5e p . Indeed we find that the scaling functions of density and 
compressibility (17.121) read off the dynamic critical exponent z = 2 and corre¬ 
lation length exponent u — 1/2 by comparing with the universal forms given 
in [5^]. In Fig. [HI we confirm these universal scaling forms though numeri¬ 
cal solution of the TBA equations (13.150 . Similarly, entropy, magnetization, 
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specific heat and susceptibility can also map out the quantum criticality of 
the model. 

8. Conclusion 

In conclusion, we have proposed a ID integrable spin-3/2 fermionic gas of 
cold atoms with spin SO( 4) symmetry. The symmetry, conserved quantities 
and integrability of this model has been studied by means of the BA We 
have shown that the integrable AO(4) symmetry spin-3/2 Fermi gas exhibits 
the spin singlet and quintet Cooper pairs in the two sets of SU (2) <g) SU (2) 
spin subspaces. The Bethe ansatz equations and finite temperature thermo¬ 
dynamical equations have been derived for the model in an analytical way. 

In particular, using the Bethe ansatz exact solutions we have thoroughly in¬ 
vestigated spin pairing, full phase diagrams, equation of states, elementary 
excitations, correlation functions, magnetism and quantum criticality of the 
model. We have also shown that the SO (4) symmetry Fermi gase possesses 
various phases of Luttinger liquids with novel a magnetism. It is particular 
interesting that in the pure paired phase, breaking a 2-string of v leads to 
four spin-1 spinons in spin excitations. The Luttinger liquid physics and 
and universal scaling behaviours of thermodynamical properties in regard 
of different spin states provide insights into understanding the large spin 
phenomena in the ID interacting fermions. 

Moreover, long distance asymptotics of various correlation functions have 
been calculated by using conformal field theory. In particular, the FFLO- 
like pair-pair correlations exist in the mixed phase of spin-2 pairs and single 
spin-3/2 atoms. Furthermore, the density profiles of the trapped gas has 
been also discussed in the context of a experimental setting with ultracold 
atoms. It turns out that the trapping center can be a mixture of different 
types of spin pairs and unpaired single atoms accompanied by the multiple 
shells involving the fully-polarized spin-2 pairs, the FFLO state, and the 
FF phase. Our results open to further study of the non -SU(k) symmetry 
interacting fermions of ultracold atoms in theory and experiment. 

Acknowledgement 

This work is supported by NSFC (Grant Nos. 91230203, 11174335,11304357, 
11434013 and 11374331), CAEP, the National Basic Research Program of 
China under Grant 2011CB922200 and 2012CB922101, and the grant from 


33 



Chinese Academy of Sciences. XWG acknowledges the Beijing Computa¬ 
tional Science Research Center and KITPC, Beijing for their kind hospitality. 
He has been partially supported by the Australian Research Council. 

References 

[1] T.-L. Ho and S. Yip, Phys. Rev. Lett. 82 (1999) 247. 

[2] S.-K. Yip and T.-L. Ho, Phys. Rev. A 59 (1999) 4653. 

[3] C. Wu, Phys. Rev. Lett. 91 (2003) 186402. 

[4] C. Honerkamp and W. Hofstetter, Phys. Rev. Lett. 92 (2004) 170403. 

[5] A. V. Gorshkov, et al, Nat. Phys. 6 (2010) 289. 

[6] M. A. Cazalilla, A. F. Ho and M. Ueda, New J. Phys. 11 (2009) 103033. 

[7] T. Lompe, et al, Science 330 (2010) 940. 

[8] J. R. Williams, et al., Phys. Rev. Lett. 103 (2009) 130404. 

[9] J. H. Huckans, et al, Phys. Rev. Lett. 102 (2009) 165302. 

[10] B. J. De Salvo, et al, Phys. Rev. Lett. 105 (2010) 030402. 

[11] X. Zhang, et al, Science 345 (2014) 1467. 

[12] S. Tale, et al, Phys. Rev. Lett. 105 (2010) 190401. 

[13] S. Tale, et al, Nat. Phys. 8 (2012) 825. 

[14] F. Scazza, et al, Nature Phys. 10 (2014) 779. 

[15] J. S. Krauser, et al, Nat. Phys. 8 (2012) 813. 

[16] M. A. Cazalilla and A. M. Rey, Progress in Physics 77 (2014) 124401. 

[17] C. N. Yang, Phys. Rev. Lett. 19 (1967) 1312. 

[18] M. Gaudin, Phys. Lett. A 24 (1967) 55. 

[19] C. Wu, Phys. Rev. Lett. 95 (2005) 266404; 

C. Wu, Mod. Phys. Lett. B 20 (2006) 1707. 


34 



[20] P. Lecheminant, E. Boulat and P. Azaria, Phys. Rev. Lett. 95 (2005) 
240402. 

[21] S. Capponi, G. Roux, P. Lecheminant, P. Azaria, E. Boulat, and S. R. 
White, Phys. Rev. A 77, 013624 (2008). 

[22] G. Roux, S. Capponi, P. Lecheminant, and P. Azaria, Eur. Phys. J. B 
68, 293 (2009). 

[23] B. Sutherland, Phys. Rev. Lett. 20 (1968) 98. 

[24] X.-W. Guan, M. T. Batchelor, C. Lee and J.-Y. Lee, Europhys. Lett. 86 
(2009) 50003. 

[25] P. Schlottmann and A. A. Zvyagin, Phys. Rev. B 85 (2012) 024535; 

P. Schlottmann and A. A. Zvyagin, Phys. Rev. B 85 (2012) 205129. 

[26] D. Controzzi and A. M. Tsvclik, Phys. Rev. Lett. 96 (2006) 097205. 

[27] Y. Jiang, J. Cao and Y. Wang, Europhys. Lett. 87 (2009) 10006; 

Y. Jiang, J. Cao and Y. Wang, J. Phys. A 44 (2011) 345001. 

[28] K. Rodrguez, A. Argiielles, M. Colome-Tatche, T. Vekua, and L. Santos, 
Phys. Rev. Lett. 105 (2010) 050402. 

[29] G. Pagano, et al, Nat. Phys. 10 (2014) 198. 

[30] S. Inouye, et a/., Nature 392 (1998) 151. 

[31] T. Bergeman, et al, Phys. Rev. Lett. 91 (2003) 163201. 

[32] M. Olshanii, Phys. Rev. Lett. 81 (1998) 938. 

[33] E. H. Lieb and W. Liniger, Phys. Rev. 130 (1963) 1605. 

[34] V. Korepin, N. M. Bogoliubove and A. G. Izergin, Quantum Inverse 
Scattering Method and Correlation Functions, Cambridge University 
Press, Cambridge, 1993. 

[35] M. J. Martins and P. B. Ramos, Nucl. Phys. B 500 (1997) 579. 

[36] P. Schlottmann, Phys. Rev. Lett. 68 (1992) 6132. 


35 



[37] M. Takahashi, Pro. Theo. Phys. 46 (1971) 1388. 

[38] C. N. Yang and C. P. Yang, J. Math. Phys. 10 (1969) 1115. 

[39] X.-W. Guan and T.-L. Ho, Phys. Rev. A 84 (2011) 023616. 

[40] S. Chen, X.-W. Guan, X. Yin, L. Guan and M. T. Batchelor, Phys. Rev. 
A 81 (2010) 031608. 

[41] L. Guan and S. Chen, Phys. Rev. Lett. 105 (2010) 175301. 

[42] Y. Liao, et al. } Nature 467 (2010) 567. 

[43] G. Orso, Phys. Rev. Lett. 98 (2007) 070402. 

[44] H. Hu, X.-J. Liu and P. D. Drummond, Phys. Rev. Lett. 98 (2007) 
070403. 

[45] X.-W. Guan, M. T. Batchelor, C. Lee and M. Bortz, Phys. Rev. B 76 
(2007) 085120. 

[46] X.-W. Guan, M. T. Batchelor and C. Lee, Rev. Mod. Phys. 85 (2013) 
1633. 

[47] J.-Y. Lee and X.-W. Guan, Nucl. Phys. B 853 (2011) 125. 

[48] F. H. L. Essler, H. Frahm, F. Gohmann, A. Klumper and V. E. Ko- 
repin, The One-Dimensional Hubbard Model, Cambridge University 
Press, Cambridge, 2005. 

[49] N. M. Bogoliubove and V. E. Korepin, Theor. Math. Phys. 82 (1990) 
231. 

[50] H. Frahm and A. Schadschneider, J. Phys. A: Math. Gen. 26 (1993) 
1463. 

[51] H. Frahm and V. E. Korepin, Phys. Rev. B 43 (1991) 5653. 

[52] X.-W. Guan, X. G. Yin, A. Foerster, M. T. Batchelor, C. H. Lee and H. 
Q. Lin, Phys. Rev. Lett. Ill (2013) 130401. 

[53] S. Sachdev, Quantum Phase Transitions, Cambridge University Press, 
Cambridge, 1999. 


36 



[54] N. Gemelke, X. Zhang, C.-L Hung and C. Chin, Nature 460 (2009) 995. 

[55] C.-L. Huang, X. Zhang, N. Gemelke and C. Chin, Phys. Rev. Lett. 104 
(2010) 160403. 

[56] C.-L. Huang, X. Zhang, N. Gemelke and C. Chin, Nature 470 (2011) 
236. 

[57] X. Zhang, C.-L. Hung, S.-K. Tung and C. Chin, Science 335 (2012) 1070. 

[58] M. P. A. Fisher, P. B. Weichman, G. Grinstein and D. F. Fisher, Phys. 
Rev. B 40 (1989) 546. 


37 



